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The calculation of the long-distance contribution to the — mass matrix is divided into three 
parts: First, the calculation of the matrix element between kaon states of the product of two 
space-time integrated, A5 = 1, four-quark weak operators. Second an RI/MOM subtraction to 
remove the short distance part of this matrix element in a fashion consistent with the calculation 
of the physical short distance part. Third an application of the Lellouch-Luscher method, gener- 
alized to second order in the weak interactions, to control finite volume errors. Such an approach 
promises to permit accurate lattice calculation of the K^-Ks mass difference and the long-distance 
contributions to Ek- 
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Norman H. Christ 



1. Introduction 

Lattice QCD has been very successful at computing the effects of the electroweak interactions 
on the properties of the strongly interacting particles. For many processes the large mass of the 
and Z bosons cause their interactions with the quarks and gluons of the hadrons to take place in a 
very small space-time region. These short distance interactions can be evaluated using electroweak 
and QCD perturbation theory and their low energy effects on hadrons described by effective four 
quark operators. For example, this approach provides a good description of both first-order decays 
and even some second order processes such as the CP violating effects in — and — B*' 
mixing. However, for general second order processes, in which two and/or Z bosons appear, 
it is possible that while each or Z exchange will appear to take place at a point, the points 
locating these two exchanges may be separated by a much larger distance ~ 1/Aqcd- Such long 
distance effects are believed to contribute to the CP violation seen in - K mixing at the 5% 
level ^ but on at least the 20% level for the CP conserving — Kg mass difference. ^ 

Here we present a method to compute such long distance effects using lattice QCD, focused on 
the case of the Ki — Ks mass difference. There are three complications which must be overcome. 
First we need to devise an Euclidean space expectation value which can be evaluated in lattice QCD 
and which contains the second order energy shift of interest. 

Second, such a lattice quantity will involve a product of two, first-order weak Hamiltonian 
densities, 2, each corresponding to one of the or Z exchanges. The short distance 

behavior of this product as |jci — .x;2 1 — will not describe the actual behavior of the exchange of 
two or Z bosons at nearby space-time points. Thus, this incorrect short distance behavior must 
be removed and replaced by the known, physical, short distance behavior described above. 

Third, the effects of finite volume, necessary in a lattice calculation, must be removed. These 
appear especially significant since the infinite volume expression contains continuous integrals, 
often with vanishing energy denominators evaluated as principal parts, while the finite volume 
quantity is a simple sum of discrete finite volume states. Here a generalization of the method of 
Lellouch and Luscher iQ] can be used. We will now discuss how each of these obstacles may be 
overcome in a calculation of the Kl — Ks mass difference, Aitik- 



2. Second order lattice amplitude 

The standard description of — Z° mixing provides an expression for the Kl — Ks mass 
difference which we will write as 

Am^ = 2^y:^^!i^^^^^M^^. (2.1) 
^ niK-Ea 

Here CP violating effects, at the 0.1% level, have been neglected, we are summing over inter- 
mediate states I a) with energy Ea and normalization factors associated with the conserved total 
momentum are suppressed. This generalized sum includes an integral over intermediate state ener- 
gies and the ^ indicates the principal part of the integral over the Ea = mn singularity. 

'For a discussion of the lattice QCD calculation of long distance effects in different decay processes see Ref. pl]. 
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One possible way to capture a similar expression in a Euclidean space lattice calculation is 
to evaluate the time-integrated second-order product that if evaluated in Minkowski space would 
yield the Am/j contribution to the time evolution over a time interval [fa,?/,]: 



dt2 / ' dtiHv,{t2)Hy^{h)K\ti)). 



(2.2) 



Here the initial is created by a source ^\ti) at the time and the final state destroyed by 
the sink ^{tf) at time tf. This amplitude is represented schematically in Fig. |l|. Equation 2.2 can 
be evaluated as a standard Euclidean space path integral with tf ^ th ^ ta ^ tf. If the time extent 
of this Euclidean path integral is sufficiently large, then when converted to an operator expression, 



Eq. becomes the vacuum expectation value of the time-ordered product of Heisenberg opera- 

— n 

tors. Assuming that tf — th and ta — tj are sufficiently large to project onto the K and A' states, 
substituting a sum over energy eigenstates \n), and integrating over t2 and ti one obtains: 



I 

n^no 

2^ 



(t\Hw\n){n\Hw\K'^) 



th — ta 



-(E„-mK)(th-ta) _ I 



-{tf-ti)inK 



mk -En rriK-En 



taf-{t'\Hw\no){no\Hw\Ky-^'f-''^'-'^. 



(2.3) 



Anticipating a result from Sec. ^ we have assumed that a single two-pion intermediate state |?io) is 
degenerate with the kaon and treated that state separately in the time integrations. 
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Figure 1: One type of diagram contributing to jz/ of Eq. 2.2. Here f2 and t\ are integrated over the interval 
[f„,ffo], represented by the shaded region between the two vertical lines. In addition to this connected quark 
flow there will also be disconnected diagrams in which no quark lines connect Hw {ti) and Hw (fi )• 



The coefficient of the {th — ta) term in Eq. 2.3 is then a finite volume approximation to Am^: 

(2.4) 



Am 



FY . 



21 



{^\Hw\n){n\Hw\K^) 



rriK - En 



The other terms in Eq. ^ fall into four categories: i) The term independent of th — ta within the 
large curly brackets. This constant must be distinguished from the desired term proportional to 
th — ta- ii) Exponentially decreasing terms coming from states \n) with > m^. These should be 
negligible if th — ta is sufficiently large, iii) Exponentially increasing terms coming from states \n) 
with En < ifiK- These will be the dominant contributions and must be accurately determined and 
removed as discussed in the paragraph below^. iv) The final term proportional to {th — taY arises 
because our choice of volume makes one 71 — 71 state, |?io)> degenerate with the kaon. 

^The author thanks Guido Martinelli and Stephen Sharpe for pointing out this behavior which had been overlooked 
when this talk was presented. 
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The exponentially growing terms pose a significant challenge. Fortunately, we have some free- 
dom to reduce their number and complexity. The two leading terms corresponding to the vacuum 
and single pion states can be computed separately and subtracted. Two pion states lying below niK 
can be eliminated using the same techniques that have been developed to evade the Maiani-Testa 
theorem and force the lowest energy k — n state to be the on-shell K ^ nn decay product. Either 
choosing the kaon to have a non-zero laboratory momentum of 753 MeV or introducing G-parity 
boundary conditions to force non-zero pion momentum can eliminate all n — n states with energy 
below niK, at least for those lattice volumes that will be accessible within the next few years. 



3. Short distance correction 

The product of operators appearing in Eq. accurately describes the second order weak ef- 
fects when the corresponding Hamiltonian densities 2 are evaluated at space-time points 
separated by a few lattice units a: \x2 — xi | ^ a. However, as |x2 — xi | ^0 the behavior is unphys- 
ical, being dominated by lattice artifacts rather than revealing the short distance structure of 
and Z exchange. Fortunately, non-perturbative Rome-Southampton methods can be applied here 
to accurately remove this incorrect behavior and replace it with the correct short distance behavior, 
that portion of the process that has been traditionally computed using lattice methods. 

This can be done by identifying the short distance part of the amplitude by evaluating the 
four-quark, off-shell Green's function 

^aPysiPi) = {da{P4)sp{p3) J d'^Xid'^X2.^w{x2)'^w{x\)Sy{p2)ds{p\)) ■ (3.1) 

Here the quark fields are Fourier transformed and the gauge is fixed. A class of connected contri- 
butions to this Green's function is shown in Fig. ^ A standard application of Weinberg's theorem 
demonstrates that if the external momenta pi obey a condition such as Pi- Pj = /x^(l —45,^), then 
for }j} ^ Aq(-.q all of the internal momenta contributing to T{pi) will have the scale jx, up to terms 
of order A^(^q//x2. 




Figure 2: Diagram representing a class of connected contributions to the off-shell, four-quark Green's 



function defined in Eq. [3.1| . In a non-perturbative evaluation of the Green's function r(p,), graphs of this 



sort including all possible gluon exchanges would be included. 

At low energies this high momentum part of the integrated product .y^(x2) J^(xi) can be 
represented as a linear combination of four-quark operators {Os}\<s<s- These operators are typi- 

in 



cally normalized by imposing conditions on off-shell Green's functions similar to that in Eq. 3.1 
which the product of .j^w operators is replaced by Os and the same kinematic point evaluated. The 
result is an alternative expression for the short distance part of the amplitude .s^ : 

^SD = ^{tf) r dxo f d\ f cf(At')0,(x,xo)/:" (?,•))• (3-2) 
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The functions c'f'(/i^) are Wilson coefficients for the lattice-regularized operator product expanded 
in operators normalized using the regularization invariant (RI) Rome-Southampton scheme. 

Thus, we can replace the incorrect short distance part of our lattice operator product by the 
correct continuum contribution by adding to the integrated operator product in Eq. the operator: 

r^^o / J^^f {cr'(M')-cr(M')}0.v(x,xo). (3.3) 

'n .V=l 

Here the {c™"'}i<i<5 are the usual continuum Wilson coefficients that are computed from elec- 
troweak and QCD perturbation theory to represent the correct short distance part of the physical 
second order weak process while the lattice coefficients c'f' can be computed from the somewhat 



elaborate but well defined lattice RI/MOM calculation of the Green's functions in Eq. |3.1. 

An important issue on which the above argument depends is the degree to which the dimension- 
6, four-quark operators introduced above capture the entire short distance part of the lattice ampli- 
tude. Since the degree of divergence of the diagram shown in Fig. ^ is +2, the integration that 



remains after the "subtraction" of the operator in Eq. ^ will still receive 0(1) contributions from 
the lattice scale. This difficulty can be avoided by including dimension eight terms in the Wilson 
expansion employed in Eq. A more physical and more practical approach includes the charm 



quark in the lattice calculation so that GIM suppression makes the integration more convergent. 



4. Controlling finite volume errors 

We now turn to the heart of this proposal: a demonstration that the potentially large volume 



dependence coming from those energy denominators in Eq. |Z3| with £"„ ~ niK can be removed, 
leaving C?(l/L^) finite- volume errors. This important conclusion is a consequence of a generahza- 
tion of the original method of Lellouch and Luscher. The starting point is Luscher's relation 
between an allowed, finite- volume, two-particle energy, E = and the two-particle scat- 

tering phase shift 5{E): 

^{kL/2n) + 5{E) = nn (4.1) 

where n is an integer and the known function <^{q) is defined in Ref. [Q]. 

Following Lellouch and Luscher we consider the i'-wave, 7i — 7i scattering phase shift as mod- 



ified by the weak interactions and use Eq. ^^Jj to connect this to the finite volume energies, deter- 
mined using degenerate perturbation theory for the Kg-^ n — n finite volume system. For simplicity 
we will Umit our discussion to the larger A/ = 1/2 part of Hw and the I = n — 7i state.^ The rela- 



tion between the finite and infinite volume second order mass shift is obtained by imposing Eq. 4-. 1 
accurate through second order in Hw- 

We begin by examining the energies, accurate through second order in the strangeness chang- 
ing, A/ = 1/2 weak Hamiltonian, Hw, of the finite volume system made up of a Ks meson, an / = 
two-pion state \no) with energy E„g nearly degenerate with niK and other single and multi-particle 
states coupled to and \no) by Hw- Following second order degenerate perturbation theory, we 



^Treating the general case is not difficult: ^^"^ -H^ ^''^ and '''^ ■ '''^ can be analyzed in the same 



way while the combination //^ '''^ -H^^^^^ contains no two-pion intermediate states. 
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can obtain the energies of the Ks and two-pion state |«o) as the eigenvalues of the 2x2 matrix: 

(4.2) 

Finite and infinite volume quantities can then be related by requiring that the eigenvalues of 



the 2x2 matrix in Eq. |4^ solve Eq. |4jJ where the phase shift 5{E) is the sum of that arising 
from the strong interactions, 8q{E), a resonant contribution from the Ks pole and more familiar 
second-order Born terms: 



= 5o(£') + arctan( 



r(£)/2 



niK + Am^ — E ' 



\{mw\no)? 
E-Ea 



(4.3) 



Here r(£') is proportional to the square of the Ks - two pion vertex which becomes the Ks width 
when evaluated at £" = niK- 



riE) = 27i\{KK{E)\Hw\Ks)\\ 



(4.4) 



where for the infinite volume, / = 0, s-wave, 2-pion state we choose the convenient normalization 
{n7r{E)\nn{E')) = 8{E — E'). The three terms in Eq. are shown in Fig. |3[ 





Figure 3: Diagrams showing the three contributions to the n — n phase shift when both strong and second 
order weak effects are included. The states /3 are multi-particle states with 5 = ±1. 

The easiest case to examine is that in which the volume is chosen to make — thk very small 
on the scale of Aqcd but large compared to Y or Am^, so that niK and En^ are not "degenerate". 



Expanding Eq. ^ and the K — n energy eigenvalue from Eq. ^ in Hyv and collecting all terms of 
second order in H\f we find: 



dE 



\{Ks\Hw\no){' 



E=E„. 



niK 



n^Ks 



{n\Hw\nQ)Y 
- F 



-no 



-no 



-niK 



-no 



-no 



This relation has two useful consequences. First we can equate the residues of the kaon poles, 
niK on the left- and right-hand sides. This gives us the original Lellouch-Luscher relation. 



Second we can subtract the pole terms and equate the remaining parts of Eq. evaluated at 
Eno = ^K- This second result will be used below to remove the second-order Born terms. 

Finally, closer to the original spirit of Lellouch and Luscher, we substitute the phase shift 8{E) 



from Eq. 4.3 into Eq. 4.1 and require that the resulting equation be valid at the energy eigenvalues 
E± of the 2x2 matrix in Eq. 4.2 for a box chosen to make = m/f. To zeroth order in H\f, 
this relation is the usual Luscher relation between 5q{E) and the allowed, finite volume, n — n 



energy. When Eq. 4.1 is expanded to first order, we reproduce the standard derivation of Lellouch 
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and Luscher's relation. Expanding to second order in Hw yields the desired relation between the 
finite and infinite volume expressions for the second order weak contribution to the Ks difference: 



Am 



Ks 



\{n\H^^\Ks) 



niK 



+ ■ 



1 



i52(0 + 5o) 

2 dE^ 



\{no\Hw\Ks 



dE, 



no 



dE 



\{n\Hw\Ks)f 



E=E, 



"0 



(4.6) 



where Eq. [4.5[ , evaluated at En^ = tuk with the pole terms subtracted has been used to eliminate the 
second-order Bom terms. Note the d /dEn^ appearing in the final term in Eq. 4.6 must be evaluated 
by varying the spatial volume which determines En^^.^ 

To obtain the Ki — Ks mass difference we first observe that the Ki second order mass shift is 
given by a formula similar to Eq. ^|6|in which Kl replaces Ks and all but the first term on the right- 
hand side is omitted since Kl does not couple to two pions, assuming CP conservation. Second if 



this new equation is subtracted from Eq. |4^ the result is similar to Eq. |4^ with Am = Amk^ — Amx^ 
on the left-hand side, the first term on the right-hand side is simply Am^ of Eq. lA and the 
remaining two 0(1 /Lr') correction terms on the right hand side of Eq. are unchanged. 



5. Conclusion 

We have proposed a lattice method to compute the Kl — Ks mass difference in which all errors 
can be controlled at the percent level. Both short and long distance effects are represented, includ- 
ing a possibly A/ = 1 /2-enhanced contribution from 7 = two pion states. Given the complexity of 
the analysis, the importance of physical kinematics and the difficulty of the disconnected diagrams 
this calculation is not practical today but may be possible in the next few years. 

The author thanks his RBC/UKQCD collaborators for important contributions to this work and 
Laurent Lellouch, Guido Martinelli and Stephen Sharpe for very helpful discussions. This work 
was supported in part by U.S. DOE grant DE-FG02-92ER40699. 
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In a one-dimensional example, these derivative terms come naturally from a generalization of the usual contour 
integration relation between finite volume sums and infinite volume momentum integrals which includes the effects of a 
double pole arising from the vanishing on-shell energy denominator. 
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